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TECHNICAL PUBLICATION 


ESTIMATING THE MEAN ANNUAL SURFACE AIR TEMPERATURE 
AT ARMAGH OBSERVATORY, NORTHERN IRELAND, AND THE GLOBAL 
LAND-OCEAN TEMPERATURE INDEX FOR SUNSPOT CYCLE 24, 
THE CURRENT ONGOING SUNSPOT CYCLE 


1. INTRODUCTION 


As noted by Gray et al., Sir William Herschel was the first to suggest a possible close connec- 
tion between the Sun and the Earth’s climate. 1-7 The Sun, being the source of energy that impacts 
and drives the Earth’s climate system, displays a variety of changes over both short and long term 
time scales, the most obvious examples being the somewhat regular waxing and waning of sunspots 
with time (i.e. , the sunspot cycle (SC)), first described by Samuel Heinrich Schwabe, 8 a German 
apothecary and amateur astronomer who observed the Sun from Dessau, Germany, 9-1 3 and the now 
well established variation of the Sun’s irradiance over the SC. 14-24 

Other factors related to the SC have been linked to changes in climate as well. Some of these 
other factors include the role of cosmic rays and the solar wind (i.e., the geomagnetic cycle) on 
climate, 25-43 as well as the apparent close association between trends in global and northern hemi- 
spheric temperature and the length of the SC, 44-52 although some investigators have described the 
inferred association between climate and, in particular, SC length as now being weak. 53-60 

More recently, Solheim et al. have reported on the relation between SC length and the average 
temperature in the same and immediately following SC for a number of meteorological stations in 
Norway and in the North Atlantic region. 61 They noted that while they found no significant trend 
(correlation) between SC length and the average temperature when measured for the same cycle, in 
contrast, they found a significant negative trend when SC length was compared with the following 
cycle’s average temperature. From this observation, they suggested that average northern hemispheric 
temperature during the present ongoing SC (SC24) will be lower by about 0.9 °C than was seen in 
SC23 (spanning 1996-2007, based on yearly averages of sunspot number (SSN), and onset for SC24 
occurring in 2008). 

The purpose of this Technical Publication (TP) is to examine the annual variations of the 
Armagh surface air temperature (ASAT) and the Global Land-Ocean Temperature Index (GLOTI) 
in relation to SSN and the SC in order to determine their likely values during SC24. Hence, it may 
provide insight as to whether solar forcing of global temperature is now lessening as a contributor to 
global warming, thereby indicating a possible cooling in the near term immediate future that poten- 
tially could ameliorate the effect of increased anthropogenic warming. 57 ’ 58 ’ 62-66 


2. RESULTS 


2.1 Annual and 10- Year Moving Averages of ASAT, GLOTI, and SSN 

Figure 1 displays the annual (January-December; thin jagged line) and 10-year moving aver- 
ages (10-yma) (thick smoothed line) of (a) ASAT, (b) GLOTI, and (c) SSN for the interval 1840-201 1 
(although annual values for these parameters are now known through 2012, it should be noted that 
at the time of writing, the data were available only through September 2012, which forms the basis 
for this study). For ASAT, the values are available for the years 1844-2011, while for GLOTI the 
values are available for the years 1880-2011. For SSN, the values are available for the entire inter- 
val 1840-2011; however, they are considered most reliable only for the interval 1882-2011. Prior to 
1882, the values of SSN are considered less reliable based on a comparison of Zurich or Wolf SSN 
(now International SSN) against group SSN. 11 ’ 67 - 73 In each subpanel, the mean, standard deviation 
(, sd ), and low- and high-annual values are given, as well as their years of occurrence. Individual SCs 
are identified across the bottom in figure 1(c). Hence, for their individual intervals of data avail- 
ability, ASAT is found to average about 9.25 ±0.55 °C (i.e., the mean± 1 sc/ interval), GLOTI about 
-0.02 ±0.27 °C, and SSN about 55 ±43.4, and the ranges of annual values are found to span 7.4 (in 
1879) to 10.6 °C (in 2007) for ASAT, -0.42 (in 1907) to 0.63 °C (in 2010) for GLOTI, and 1.4 (in 1913, 
minimum amplitude (Rmin) for SC15) to 190.2 (in 1957, maximum amplitude (Rmax) for SC19) for 
SSN. The peak annual values of ASAT and GLOTI are about 2.4 higher than their respective long 
term means, while the peak annual mean for SSN is about 3 sd higher than its long term mean. The 
lows are respectively about 3.4, 1.5, and 1.2 Slower than their long term means. 

The ASAT is one of the longest, continuous, thermometer-based temperature records avail- 
able for study. 45 ’ 74-79 The Armagh Observatory lies about 1 km northeast of the center of the ancient 
city of Armagh, Northern Ireland, being located at 54°21'12" N. and 6°38'54" W. It is situated about 
64 m above sea level at the top of a small hill in an estate of natural woodland and parkland that 
measures about 7 ha. Previous studies have shown that its rural environment has ensured that the 
temperature measurements suffer little or no urban microclimatic effects and that the measured tem- 
perature record can serve as a good proxy for monitoring long term trends in both northern hemi- 
spheric and global annual mean temperature. 76 ’ 80 The data for the ASAT are available online at 
<http : //climate, ar m. ac . uk/scans > . 

The GLOTI is a measure of the anomaly in global land-ocean temperatures relative to the 
base period of 1951-1980, where the data are taken from the Global Historical Climate Network, 
version 3, using elimination of outliers and homogeneity adjustment. The data (January-December 
averages) for the GLOTI are available online at <http://data.giss.nasa.gov/gistemp>. 81 ” 83 
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Figure 1. Annual (thin jagged line) and 10-yma (thick smoothed line) of 

(a) ASAT for the interval 1844-2011, (b) GLOTI for the interval 
1880-2011, and (c) SSN for the interval 1840-2011. 


The SSN is a measure of the strength of the solar activity, which varies over a period (or 
SC length) of about 10-12 years, based on annual or yearly averages of SSN during the modern 
era of sunspot observations (SC 10-present). The SSN values shown in figure 1(c) correspond to 
the last 3 years of SC8 through the first 4 years of SC24, and the SSN data are available online at 
<ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUNSPOT_NUMBERS/INTERNATIONAL>. 

The 10-yma values are shown to indicate trends in the annual values. For AS AT, its 10-yma 
value was lowest in 1883 (having a value of 8.44 °C) and highest in 2002 and 2003 (having a value of 
10.13 °C). Hence, in terms of its 10-yma values, the ASAT has increased 1.69 °C in about 120 years, 
suggesting a rise in temperature at Armagh Observatory, Northern Ireland, at the average rate of 
about 0.014 °C per year or about 0.14 °C per decade (equivalent to a rise of about 0.025 °F per year or 
0.25 °F per decade). It is interesting that prior to 1883, the 10-yma value at Armagh is found to have 
been trending downwards from a high of 9.46 °C in 1850 to the low of 8.44 °C in 1883, a decline of 
about 0.031 °C per year, or 0.31 °C per decade. Also interesting is that the rise in temperature from its 
low in 1883 rose quickly to 9.1 °C in 1901, leveled off for the next 35 years, rose once again to about 
9.54 °C in 1945, and then fell to a local low of 9.05 °C in 1982 before rising once again to its highest 
observed peak temperature of 10.13 °C in 2002 and 20 03. 84 Values since 2003 have fallen slightly to 
10 °C (in 201 1). The trend line, as determined using the 10-yma values, was below its long term mean 
between about 1857 and 1932 and again between about 1964 and 1984, while it was above its long term 
mean prior to 1857, between 1933 and 1963, and after 1984 (to the present). 

For GLOTI, its 10-yma value was lowest in 1907 and 1908, measuring -0.34 °C, and highest 
in 2006, having a value of 0.55 °C. Hence, the GLOTI has increased 0.89 °C over about 103 years, 
indicating an average rise in warming of about 0.009 °C per year or 0.09 °C per decade (equivalent to 
an average rise of 0.016 °F per year or 0.16 °F per decade). Like the ASAT, it too shows a local peak 
(and flattening) in the early 1940s, followed by a rise after about 1967 to its observed peak in 2006, 
rising much faster over the past 40 years at an average rate of about 0.015 °C per year or 0. 15 °C per 
decade. The trend line in the GLOTI was always below its long term mean prior to about 1940, while 
it was near its long term mean from about 1940-1970, and it has always been above its long term mean 
thereafter. 

For SSN, its 10-yma values appear to have been trending downwards between about 1850 and 
1900, rising slightly but remaining below its long term mean until about 1935, and then rapidly rising 
to its peak in the late 1950s (associated with SC 19), staying above its long term mean from about 1935 
until 2003. 85-88 The trend line is now below its long term mean and is expected to stay below it for the 
near term foreseeable future. While the trend line is currently below its long term mean, the annual 
SSN value rose above it in 201 1 and will remain above it briefly, probably through about 2013-2015 
since SC24 is now in its period of maximum phase, which typically lasts about 3-4 years in length. 

Figure 2 shows the residuals for (a) ASAT and (b) GLOTI, where the residual is simply the 
difference between the observed annual value and the corresponding 10-yma value (i.e., for the same 
year). For ASAT, the mean is 0 °C and the sd is 0.41 °C. Hence, use of the 10-yma values reduces 
the variance by about 45%. Furthermore, runs-testing of the residuals yields a normal deviate of the 
sample z = 0.96, which by hypothesis testing suggests that the residuals are randomly distributed. 89 
Therefore, use of the 10-yma values provides a convenient way to determine the overall yearly trend in 
the annual values of the ASAT. 
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Figure 2. Residuals for (a) ASAT and (b) GLOTI, where the residual is 

simply the difference between the observed annual values and the 
corresponding 1 0-yma values (for the same year). 
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Similarly, for GLOTI, the mean is 0 °C, and the sd is 0.08 °C, yielding a reduction of the vari- 
ance by about 90%. Runs-testing of its residuals also are found to be randomly distributed (z = -1 .08). 
Hence, as with the AS AT, the 10-yma values provide a convenient way to determine the overall yearly 
trend in the annual values of the GLOTI. 

Figure 3 depicts the distributions of the residuals for (a) ASAT and (b) GLOTI. About 
80% of the time, the 10-yma values for ASAT and GLOTI are within the interval 0±0.5 °C and 
0 ± 0. 1 °C, respectively. Hence, given the annual values of ASAT and GLOTI, one can infer approxi- 
mately their yearly 10-yma values (up to 5 years beyond the last available 10-yma value). For the year 
2007 then, one expects the 10-yma value for ASAT to be about 10.6 ± 0.5 °C and GLOTI to be about 
0.59 ±0.1 °C, inferring only about a 10% chance that the 10-yma values will fall below 10.1 °C and 
0.49 °C, respectively. 



Figure 3. Distributions of the residuals for (a) ASAT and (b) GLOTI. 


Figure 4(a) shows the scatter plot of the annual values of GLOTI versus ASAT for the interval 
1880-201 1. The thin vertical and horizontal lines represent the median values of the parameters, and 
the thick diagonal line represents the inferred preferential regression of GLOTI versus ASAT. Hence, 
as the annual value of the ASAT increases (decreases), one expects the annual value of the GLOTI 
to increase (decrease). The inferred regression is given by y = -2.822 + 0.302x, having a coefficient 
of linear correlation (r) = 0.6 and coefficient of determination (r 2 ) = 0.36, meaning that about 36% 
of the variance in the GLOTI annual values can be explained by the yearly variation in the ASAT 
annual values (or vice versa) during the interval 1880-2011. The inferred regression has a standard 
error of estimate ( se ) of 0.22 °C, and hypothesis testing of the slope reveals that the inferred regres- 
sion is statistically important at confidence level (c/)> 99.9%. Also shown in figure 4(a) is the result 
of Fisher’s exact test for the observed 2x2 contingency table (and for any other result that shows 
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Figure 4. Scatter plots of (a) annual values of GLOTI versus AS AT for the interval 1880-201 1 
and (b) GLOTI versus AS AT using 10-yma values for the interval 1885-2006. 


a greater departure from independence 90 ) being the probability (P ) « 0. 1%. Hence, there is strong sta- 
tistical evidence linking ASAT and GLOTI, with about 70% (= 93/132) of the annual values obeying 
the paradigm of warmer (cooler) values of GLOTI being associated with warmer (cooler) values of 
ASAT 


However, a peculiarity is noticeable in figure 4(a); in particular, the annual values for the 
year 2010. For GLOTI, the annual value for 2010 measures 0.63 °C, the warmest index value on 
record; while for ASAT, the annual value for 2010 measures only 8.75 °C, the coolest temperature 
measured at Armagh since 1986. For the first four months of the year 2010, the Oceanic Nino Index 
(ONI) indicated the presence of a weakening El Nino (EN), followed by two months of neutral 
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conditions and then a strengthening La Nina (LN) for the remainder of the year. In contrast, ONI 
for the year 1986 displayed neutral conditions for the hrst seven months of the year followed by 
a strengthening EN for the remainder of the year (ONI monthly anomalies are available online at 
<http://www.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml>). So, it is 
somewhat puzzling as to why such a large discrepancy exists between the inferred GLOTI and the 
observed GLOTI for the year 2010 given the observed AS AT. The discrepancy measures 3.7 se based 
on the inferred, highly statistically important positive preferential regression determined for GLOTI 
versus ASAT. 

Figure 4(b) displays the scatter plot of GLOTI versus ASAT using 10-yma values. The scatter 
plot strangely suggests two different regimes of regression hts: one fit active for the interval 1885— 
1945 and the other fit active from about 1982-2006 (the year of the last available 10-yma values used 
in this TP), with a transition occurring between them. Both regressions are inferred to be statisti- 
cally important, with the regression for the current interval having r = 0.962, r 2 = 0.925, se = 0.036 °C, 
and cl> 99.9%. While true, the observed 10-yma value of GLOTI for 2006 (0.55 °C) is about 0.1 °C 
warmer than that found using the inferred regression for 1982-2006, and it is about 0.4 °C warmer 
than that found using the 1885-1945 regression given the 10-yma value of ASAT for 2006 (10 °C). 
One can only speculate as to why there exists this observed discrepancy represented by two statisti- 
cally important inferred regressions with a transition occurring between them. For example, are the 
discrepancies due to the increased concentration of anthropogenic gases in the Earth’s atmosphere, 
thereby altering the relationship between GLOTI and ASAT; or are the discrepancies somehow 
related to the changing phase of the Atlantic Multidecadal Oscillation (AMO); or are the discrepan- 
cies due to some unknown problems with one or both of the indices; or are the discrepancies due to 
some unknown combination of nonlinear effects; or etc.? 

Figure 5 shows the variation in first difference (fd) values of the 10-yma values for (a) ASAT 
and (b) GLOTI, where the fd is defined as the difference between the 10-yma value for year (t) + 1 
minus the 10-yma value for year t, and figure 6 shows the distribution of the fd values for (a) ASAT 
and (b) GLOTI. For ASAT, its fd values have a mean = 0 °C and sd= 0.05 °C, with the interval fd - 0 
± 0.05 °C capturing about 73% of the fd values. For GLOTI, its fd values have a mean = 0.01 °C and 
5fi/=0.01 °C, with the interval /c/=0.01 ±0.01 °C capturing about 74% of the fd values. Hence, using 
the mean and sd in the fd values allows one to approximate the next yearly 10-yma value (i.e., for 
2007). 
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Figure 5. Variation in fd values of the 10-yma values for (a) ASAT and (b) GLOTI. 
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Figure 6. Distribution of the fd values for (a) ASAT and (b) GLOTI. 
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For AS AT, its 10-yma value for 2006 measured 10 °C, inferring that its 10-yma value for 2007 
should be about 10 + 0 ±0.05 °C, or 10 ±0.05 °C, presuming that the fd value for 2007 falls within 
the mean± 1 ^prediction interval (a value slightly lower than that deduced using the yearly residu- 
als (see fig. 3 discussion). Using this value, one estimates the annual AS AT for the year 2012 to be 
about 9.58 ± 1 °C, or below 10.58 °C. It measured 10.28 °C in 2011, and the January-September 2012 
ASAT average is about 10.7 °C. (Now that the year 2012 is complete, ASAT for the year is known to 
be 9.7 °C.) 

For GLOTI, its 10-yma value for 2006 measured 0.55 °C, inferring that its 10-yma value 
for 2007 should be about 0.55 + 0.01 ±0.01 °C or 0.56 ±0.01 °C. Using this value, one estimates the 
annual value of GLOTI for the year 2012 to be about 0.69 ±0.2 °C, or below 0.89 °C. The warm- 
est GLOTI to date is 0.63 °C in 2010. The GLOTI measured 0.51 °C in 2011, and the January- 
September 2012 GLOTI average is about 0.5 °C. Because ONI presently is showing signs of a slight 
warming (being 0.6 °C in October), perhaps suggesting the possible start of an EN event in late 2012, 
one expects the GLOTI for the year 2012 to be >0.5 °C. (Again, now that the year 2012 is complete, 
GLOTI for the year is known to be 0.56 °C.) 

2.2 ASC, PER, Rmax, <SSN>, <ASAT> and <GLOTI> for SC10-SC23 

Figure 7 shows the variation of both the actual cyclic mean value (the thin dotted line) and 
the two-cycle moving average (the thick smoothed line for SC10-SC23 of (a) ascent duration (ASC) 
(i.e., the ascent period, defined herein as the elapsed time from sunspot Rmin year to sunspot Rmax 
year in the same cycle, based on annual averages of SSN), (b) the period (PER) (i.e., the SC length, 
or period, defined herein as the elapsed time between two consecutive sunspot minima years based 
on annual averages of SSN), (c) Rmax (i.e., the maximum amplitude, or maximum yearly value of 
SSN, for an SC), (d) <SSN> (i.e., the average SSN over the SC, from sunspot minimum year to the 
following sunspot minimum year), (e) <ASAT> (i.e., the average ASAT over the SC, from sunspot 
minimum year to the following sunspot minimum year), and (f) <GLOTI> (i.e., the average GLOTI 
over the SC from sunspot minimum year to the following sunspot minimum year). The use of the 
symbols <> simply indicates that the computed values are averages of the yearly parametric values 
over each SC (from sunspot minimum year to sunspot minimum year). 
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Figure 7. Variation of both the actual cyclic value (the thin dotted line) 
and the two-cycle moving average (the thick smoothed line) 
for SC10-SC23 of (a) ASC, (b) PER, (c) Rmax, (d) <SSN>, 
(e) <ASAT>, and (f) <GLOTI>. 
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On the basis of yearly SSN, the modern era of SCs is found to have ASC periods averaging 
about 3.8 ±0.7 years (i.e., the mean± 1 sd interval), ranging incrementally 3 to 5 years in length. 
Cycles with ASC = 3 years are considered fast risers, while those with ASC = 4 years are considered 
average risers, and those with ASC = 5 years are considered slow risers. For the modern era, there 
have been five SCs classified as fast risers, including SCI 1, SC 18, SC 19, SC21, and SC22; seven cycles 
classified as average risers, including SC 10, SC13-SC15, SC 17, SC20, and SC23; and only two cycles 
classified as slow risers, including SC12 and SC16. SC24, the current ongoing SC, had its Rmin in 
terms of annual SSN in the year 2008 and is presently in its fifth year of rise (2013). Presuming SC24 
is not a statistical outlier, one expects that its maximum annual amplitude will occur in either 2012 
or 2013 and that SC24 will be classified as a slow riser. (SSN measured 57.5 in 2012. Through May 
2013, SSN averages 62.) 

Concerning PER, its mean+ 1 interval is computed to be 10.9 ±0.9 years based on annual 
averages of SSN, ranging incrementally 10 to 12 years. Cycles with PER= 10 years are considered 
short-period cycles (they include SC 15, SC 16, SC 18, SC 19, SC21, and SC22), those having PER 
= 11 years are considered average-period cycles (they include SC10-SC12, and SC17), and those hav- 
ing PER= 12 years are considered long-period cycles (they include SC13, SC14, SC20, and SC23). 
Presuming SC24 is not a statistical outlier, one expects SC24 to end before the year 2020 (i.e., in 2017 
if it is a short-period cycle, in 2018 if it is an average-period cycle, or in 2019 if it is a long-period 
cycle). During the modern era, there has never been an SC having PER >12 years based on annual 
averages of SSN; so, if SC24 continues through 2020 or later, it would be a statistical outlier with 
respect to the 14 cycles comprising the modern era of SCs. (Based on the 12-month moving average 
of monthly mean SSN (also called smoothed monthly mean SSN), SC lengths are found to span 
116 months for SC22 to 151 months for SC23.) 91-93 

Concerning Rmax, its mean+ 1 sd interval is computed to be 116 ±38.4, based on annual 
averages of SSN. Cycles having Rmax >1 16 include SCI 1 and SC18-SC23, while cycles having Rmax 
<116 include SC10-SC17 and SC20. If, instead of using the mean Rmax, one chose to sort the cycles 
using the median Rmax value, then SC 17 would be grouped with the larger Rmax group. The largest 
Rmax is 190.2 (measured for SC19), and the smallest Rmax is 63.5 (measured for SC14; SC12 had 
Rmax = 63.7). SC24 is expected to have Rmax <116, probably near Rmax = 70. 94-97 

Concerning <SSN>, its mean+ 1 sd interval is computed to be 56.4 ± 19.6 based on annual 
averages of SSN for each SC. Cycles having <SSN> >56.4 include SC18-SN23, while cycles having 
<SSN> <56.4 include SC10-SC17. Again, if instead of using the mean <SSN> as the discrimi- 
nator for groupings of <SSN>, one chose to use the median <SSN> value, then SC17 would be 
grouped with the larger <SSN> group. Because of the strong behavioral similarity between <SSN> 
and Rmax (see below), one expects <SSN> to be <56.4 for SC24. 

Concerning <ASAT>, its mean + 1 ^interval is computed to be 9.23 ± 0.37 °C. Cycles having 
<ASAT> >9.23 °C include SC10 and SC17-SC23, while cycles having <ASAT> <9.23 °C include 
SC11-SC16 and SC21. SC23 (spanning 1996-2007), having <ASAT>= 10.08 °C, has the warmest 
<ASAT> of any modern era SC, being about 2.3 sd warmer than the mean. In contrast, SC 12 (span- 
ning 1878-1888), having <ASAT> = 8.46 °C, has the coolest <ASAT> on record, being about 2. 1 sd 
cooler than the mean. Based on the average observed change in <ASAT> between two consecutive 
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cycles, one estimates that <ASAT> for SC24 likely will be about 10.08 + 0.06 ±0.33 °C, or about 
10. 14 ± 0.33 °C (i.e., <ASAT> is expected to be <10.47 °C for SC24; any <ASAT> value warmer 
than 10.08 °C would be a new record high value). 

Concerning <GLOTI>, its mean+ 1 sd interval is computed to be -0.04 ±0.25 °C. Cycles 
having <GLOTI> >-0.04 °C (i.e., warmer than -0.04 °C) include SC17-SC23, while cycles hav- 
ing <GLOTI> <-0.04 °C (i.e., cooler than -0.04 °C) include SC12-SC16. SC23, having <GLOTI> 
=0.49 °C, has the warmest <GLOTI> on record, being about 2.1 sd warmer than the mean. In 
contrast, SC14 (spanning 1901-1912), having <GLOTI> =-0.32 °C, has the coolest <GLOTI> on 
record, being about 1.1 sd cooler than the mean. Based on the average observed change in <GLOTI> 
between two consecutive cycles, one estimates that <GLOTI> for SC24 likely will be about 0.49 + 0.07 
±0.09 °C, or about 0.56 ±0.09 °C (i.e., <GLOTI> <0.65 °C for SC24, and any <GLOTI> value 
warmer than 0.49 °C would be a new record high value). 

Based on the two-cycle moving average (shown to reflect the variation as related to the Hale 
cycle 98 ), it seems very likely that ASC, PER, <ASAT>, and <GLOTI> will trend upwards for 
SC23, while Rmax and <SSN> will trend downwards. In other words, it now seems quite clear that 
<ASAT> and <GLOTI> will no longer consistently track the trend in solar activity (i.e., Rmax or 
<SSN>), as they appear to have done prior to about SC21. (There has been consistent warming 
since SC21, although solar activity has been in decline, thus indicating a lessening in solar forcing on 
climate.) 99-101 

Regarding the grouping of SCs as fast risers, average risers, slow risers, short-period cycles, 
average-period cycles, and long-period cycles, one computes the following means and sds: 

(1) Fast risers (ASC = 3 years) have PER = 10.2 ± 0.45 years (the mean ± 1 sd prediction inter- 
val), Rmax= 158.8 ±19, <SSN> = 75.9 ± 17.5, <ASAT> = 9.32 ±0.21 °C, and <GLOTI> = 0.09 
±0.13 °C. 

(2) Average risers (ASC = 4 years) have PER= 1 1.43 ±0.79 years, Rmax = 98.3± 19.1, 
<SSN> = 47.7 ± 10.2, <ASAT> = 9.3 ±0.38 °C, and <GLOTI> = -0.07 ±0.3 °C. 

(3) Slow risers (ASC = 5 years) have PER= 10.5 ±0.71 years, Rmax = 70.8 ± 10, <SSN> = 37.8 
±4.5, <ASAT> = 8.76 ±0.42 °C, and <GLOTI> = -0.24 ±0.11 °C. 

(4) Short-period cycles have ASC = 3.5 ± 0.84 years, Rmax = 1 39.4 ± 40.9, <SSN> = 69.6 ±21 .9, 
<ASAT> = 9.29 ±0.21 °C, and <GLOTI> = -0.01 ±0.19 °C. 

(5) Average-period cycles have ASC = 4 ± 0.82 years, Rmax = 103.2 ± 3 1.7, <SSN> = 46. 7± 8.7, 
<ASAT> = 9.05 ±0.42 °C, and <GLOTI> = -0.17 ±0.2 °C. 

(6) Long-period cycles have ASC = 4±0 years, Rmax = 93.5 ±24.5, <SSN> = 46.3± 13.5, 
<ASAT> = 9.33 ± 0.52 °C, and <GLOTI> = -0.03 ± 0.37 °C. 
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Hence, fast risers tend to have higher Rmax, <SSN>, and <GLOTI> as compared to aver- 
age or slow risers, and short-period cycles tend to have higher Rmax and <SSN> as compared to 
average- or long-period cycles. Higher or lower values of <ASAT> do not appear to differentiate 
between ASC or PER classes (in particular, fast versus slow riser or short- versus long-period cycles), 
based on the t statistic for independent samples, used for comparing two means from independent 
samples. 89 Likewise, higher or lower values of <GLOTI> do not appear to differentiate between 
PER classes of SCs. (Presently, SC24 is believed to be a slow riser of possibly long-period; hence, one 
expects it to have Rmax and <SSN> below respective long term means, peaking in 2012 or 2013 and 
ending in 2019 or later, in terms of annual averages of SSN.) 

In their analysis, Solheim et al. 61 used SC lengths determined differently than that used here, 
which is one based on the use of annual yearly averages of SSN. Depending upon the parameter used 
to determine SC length, whether it is the number of sunspot groups, the number of spotless days, 
sunspot area, the ratio of new cycle to old cycle sunspots, curve htting, etc., any number of different 
SC lengths can be determined, and those lengths can vary by several months or more. 102 ’ 103 Using 
annual averages simplifies the issue of determining SC length (i.e. , based on annual averages of SSN, 
a cycle length is determined to be of integer length, 10 years, 1 1 years, 12 years, etc.). 

As previously noted, Solheim et al. reported no strong association between the length of 
the SC and the temperature over that same SC (as also can be gleaned above using SC lengths 
determined from annual averages of SSN (i.e., the mean temperature and associated sd for short-, 
average-, and long-period SC classes)). 61 Additionally, they reported that SCs of shorter (longer) 
length were strongly associated with warmer (cooler) temperatures during the following SC. Based 
on SC lengths determined using annual averages of SSN, one hnds a similar result, although the 
result is found not to be statistically important based on the t statistic for independent samples, nei- 
ther for <ASAT> nor <GLOTI>. For example, short-period cycles (PER = 10 years) are found to 
have <ASAT> = 9.46 ± 0.34 °C and <GLOTI> = 0.1 1 ± 0.22 °C during the following SC, while aver- 
age-period cycles (PER = 1 1 years) have <ASAT> = 9 ± 0.43 °C and <GLOTI> =-0.2 ±0.16 °C, and 
long-period cycles (PER = 12 years) have <ASAT> = 9.08 ±0.08 °C and <GLOTI> = -0.16 ±0.23 °C 
during the following SC. So, while shorter (longer) period cycles appear to be followed by warmer 
(cooler) temperature, based on the means and sd of the temperatures, the differences do not appear 
to be statistically meaningful (t= 1.8 for <ASAT> and t- 1.7 for <GLOTI> based on the six short- 
period cycles and three long-period cycles of the modern era, inferring cl < 90%). 

Figure 8 shows the scatter plots of the cyclic averages of (a) <ASAT> versus Rmax, 
(b) <GLOTI> versus Rmax, (c) <ASAT> versus <SSN>, and (d) <GLOTI> versus <SSN>, where 
the numbered hlled circles refer to the SC number. In each subpanel, the vertical and horizontal 
lines indicate the median values, and the diagonal line, inferring a direct preferential correlation (i.e., 
an increase or decrease in y being associated with an increase or decrease in x), indicates the inferred 
regression ht between the y and x parameters. Also given in each subpanel are the inferred regression 
equations y, r, r 2 , se, cl, and P, where y refers to the inferred regression equation, r to the coefficient 
of linear correlation, r 2 to the coefficient of determination (a measure of the amount of variance 
explained by the inferred regression), se to the standard error of estimate, cl to the confidence level, 
and P to the probability of occurrence by chance based on the Fisher’s exact test for 2x2 contin- 
gency tables. 
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(b) Rmax (d) <SSN> 

Figure 8. Scatter plots of the cyclic averages of (a) <ASAT> versus Rmax, (b) <GLOTI> 
versus Rmax, (c) <ASAT> versus <SSN>, and (d) <GLOTI> versus <SSN>, 
where the numbered hlled-circles refer to the SC, 


Concerning the inferred correlation between <ASAT> and Rmax, the regression is found to 
be only of marginal statistical importance (cl> 90%), having r = 0.53. Based on Fisher’s exact test for 
2x2 contingency tables, the probability of obtaining the observed result, or one more suggestive of 
a departure from independence, is computed to be P = 14.3%. Therefore, the association of <ASAT> 
against Rmax must be considered weak. Since Rmax for SC24 is anticipated to be less than the 
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median value, there appears to be only about a 30% chance that SC24 will have < ASAT> larger than 
its median based on its expected Rmax alone (from the 2x2 contingency table). 

Concerning the inferred correlation between <GLOTI> and Rmax, the regression is found 
to be statistically important (cl >95%), having r = 0.61. Also, based on Fisher’s exact test for 2x2 
contingency tables, the association of <GLOTI> against Rmax is found to be considerably stron- 
ger, having P = 0.8% (i.e., the probability of obtaining the observed result or one more suggestive of 
a departure from independence is < 1%). Hence, given that SC24 is expected to have an Rmax less 
than the median value, one expects <GLOTI> for SC24 very probably to be less than its median 
value, again based on Rmax alone. 

Concerning the inferred correlation between <ASAT> and <SSN>, the regression is found 
to be only of marginal statistical importance (cl> 90%), having r = 0.51. However, based on Fisher’s 
exact test for 2x2 contingency tables, the probability of obtaining the observed result, or one more 
suggestive of a departure from independence, is computed to be P= 1.5%. Hence, presuming that 
SC24 will have <SSN> less than its median value, one expects <ASAT> for SC24 likewise to be less 
than its median value based on <SSN> alone. 

Lastly, concerning the inferred correlation between <GLOTI> and <SSN>, the regression is 
found to be statistically important (cl >95%), having r = 0.61. Also, based on Fisher’s exact test for 
2x2 contingency tables, the association of <GLOTI> against <SSN> is found to be considerably 
stronger, having P = 0.8%. Hence, given that SC24 is expected to have <SSN> less than the median 
value, one expects <GLOTI> for SC24 very likely to be less than its median value based on <SSN> 
alone. 


However, the locations of SC 12 and SC23 (the coolest and warmest <ASAT> values, 
respectively) in the scatter plots involving <ASAT> are troubling. Based on figure 8(a), SC 12 has 
an <ASAT> value about 1.5 sc below the inferred regression value, while SC23 has an <ASAT> 
value about 2.5 se above the inferred regression value. Based on figure 8(b), SC 12 has an <ASAT> 
value about 1.8 se below the inferred regression value, while SC23 has an <ASAT> value about 2.5 se 
above the inferred regression value. 

Likewise, while SC 12 has <GLOTI> values close to its inferred values (given Rmax and 
<SSN>), the <GLOTI> values for SC23 are considerably above them (given its Rmax and <SSN>), 
being about 2.4 se above its inferred regression value based on Rmax and 2.7 se above its inferred 
regression value based on <SSN>. Therefore, something other than the Rmax or average SSN over 
the cycle appears to be driving the increased temperature anomaly as measured using <GLOTI>. 

Figure 9 displays the scatter plots of (a) <GLOTI> versus <ASAT> and (b) <SSN> versus 
Rmax. Both scatter plots are found to be highly statistically important (cl >99. 9%), with the former 
relation having r = 0.86 and the latter relation having r = 0.94. Based on Fisher’s exact test, P = 0.8% 
for the former relation and P = 1 . 5% for the latter relation. Hence, when < AS AT> is above (below) its 
median value, one strongly expects <GLOTI> to be above (below) its median value. Similarly, when 
Rmax is above (below) its median value, one clearly expects <SSN> to be above (below) its median 
value. 
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Figure 9. Scatter plots of (a) <GLOTI> versus <ASAT> and (b) <SSN> versus Rmax. 


Figure 10 shows the variation of (a) SSN, (b) ASAT, and (c) GLOTI from January 1 996— 
September 2012, as well as the variation of (d) <SSN>, (e) <ASAT>, and (f) <GLOTI> for elapsed 
time in years past the sunspot minimum year for SC23 (the unfilled histograms) and SC24 (the filled 
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Figure 10. Variation of (a) SSN, (b) AS AT, and (c) GLOTI from January 1996-September 2012 
and the variation of (d) <SSN>, (e) <ASAT>, and (f) <GLOTI> for elapsed time in 
years past the sunspot minimum year for SC23 (the unfilled histograms) and SC24 (the 
filled histograms). 


histograms). Based on yearly averages of SSN, SC23’s sunspot minimum year was 1996, and SC24’s 
sunspot minimum year was 2008. Also shown are the occurrences of EN (warm anomalies) and LN 
(cool anomalies) events as determined using monthly values of ONI, and the occurrence of signifi- 
cant volcanic eruptions (the small, filled triangles), those having a Volcanic Explosivity Index (VEI) 
of 4 or higher. (Occurrences of EN and LN events and significant volcanic eruptions have sometimes 
been associated with brief excursions in northern hemispheric and global temperatures.) 104-107 
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It is interesting that EN events often appear to be associated with increases (warming) in the 
GLOTI (e.g., 2009), while LN events appear to be associated with decreases (cooling) in the GLOTI 
(e.g., 1998-2000). While conditions in the Nino 3.4 region presently remain slightly above average in 
temperature (ONI = 0.6 °C in October 2012; see <http://cpc.ncep.noaa.gov>), strictly speaking, they 
are reflective of the ENSO-neutral phase. Hence, while it appears unlikely that a fully-coupled EN 
will develop during the latter part of 2012 or in early 2013, the increased warming in the Nino 3.4 
region is suggestive that the GLOTI will remain elevated for 2012 (see <http://www.elno.noaa.gov>, 
December 2012). 

Also interesting is that every year of the first 4 years of SC24 has had higher GLOTI val- 
ues than for the same elapsed time during SC23, although <SSN> and <ASAT> have generally 
been well below that occurring in SC23 for the same elapsed time from sunspot minimum. The first 
4 years of SC24 have had an average GLOTI of about 0.54 °C, as compared to 0.24 °C for the 
first 4 years of SC23. So far this year (i.e. , for January-September 2012), GLOTI is averaging 
0.5 °C, as compared to the first nine months of 1990 of SC23 when it averaged only 0.37 °C. So, 
unless the October-December 2012 GLOTI values are extremely low (<-0.06 °C), elapsed time year 
four of SC24 also wifi have a GLOTI value higher than was seen in elapsed time year four of SC23 
(=0.37 °C). Hence, presently it appears that SC24 wifi have <GLOTI> values higher than was 
seen for SC23 (=0.49 °C), unless, of course, the latter years of SC24 have cooler GLOTI values 
(<0.46 °C, assuming SC24 has PER =12 years). (The year 2012 is now known to have AS AT =9.7 °C 
and GLOTI=0.56 °C. The ASAT value is lower in comparison to the same elapsed time for SC23, 
9.93 °C; the GLOTI value is higher in comparison to the same elapsed time for SC23, 0.4 °C.) 

Concerning the two significant volcanic eruptions occurring between 1996 and 2012, although 
both had VEI = 4, neither eruption appears to have had any major, long-lasting effect on either 
ASAT or GLOTI. The two eruptions occurred on August 15, 1999 (Shiveluch, Kamchatka, Russia, 
56.653° N., 161.36° E.) and June 4, 201 1 (Puyehue, Chile, 40.59° S., 72. 117° W.). The volcanic erup- 
tion data are available online at <http://www.ngdc.noaa.gov>. 

2.3 Behavioral Aspects of ASAT and GLOTI During Rmin and Rmax Years 

Ligure 1 1 displays the yearly averages of (a) Rmin (filled circles) and Rmax (filled triangles), 
(b) ASAT for Rmin years (filled circles) and GLOTI for Rmin years (filled triangles), and (c) ASAT 
for Rmax years (filled circles) and GLOTI for Rmax years (filled triangles) for SC9-SC24. Visu- 
ally, the behaviors of the parameters are strikingly similar, especially the one between Rmax and 
Rmin. 108 ’ 109 While solar activity now appears to be tracking downwards, in contrast, ASAT and 
GLOTI are at or near record highs for SC23 and SC24. 
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(c) Sunspot Cycle 

Figure 11. Yearly averages of (a) Rmin (filled circles) and Rmax (filled triangles), (b) ASAT for 
Rmin years (filled circles) and Rmax years (filled triangles), and (c) GLOTI for Rmin 
years (filled circles) and Rmax years (filled triangles) for SC9-SC24. 


Figure 12 shows the scatter plots of (a) Rmax versus Rmin, (b) ASAT (Rmax year) versus 
ASAT (Rmin year), and (c) GLOTI (Rmax year) versus GLOTI (Rmin year), with the arrows located 
along the v-axes marking the values in the Rmin year for SC24. Concerning figure 12(a), it is appar- 
ent that one can easily predict the size of the ongoing cycle 2-4 years in advance from the Rmin value 
alone. Using SC9-SC23, one finds the inferred regression to be statistically important at cl >95%, 
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(c) GLOTI (Rmin Year) 

Figure 12. Scatter plots of (a) Rmax versus Rmin, (b) AS AT (Rmax year) versus ASAT 

(Rmin year), and (c) GLOTI (Rmax year) versus GLOTI (Rmin year), with the 
arrows located along the x-axes marking the values in the Rmin year for SC24. 
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having r = 0.53. Applying the Rmin value for SC24 (=2.9), one deduces the ±1 se prediction interval 
for SC24’s Rmax to be 93.9 ± 32.6, inferring that SC24 probably will have Rmax >61.3, in close agree- 
ment to the predictions of de Jager and Duhau 95 and Wilson. 97 However, it is interesting that if one 
ignores SC 19 (the largest Rmax SC), the resultant inferred regression becomes even stronger, having 
y = 61.551 +6.826x, r = 0.787, r 2 = 0.619, se = 2Q.l, and cl> 99.9%, thus, yielding Rmax = 81.3±20.7 
for SC24 (i.e., Rmax >60.6). As previously noted, SSN for 201 1 measures 55.7; hence, one anticipates 
2012 to be the Rmax year for SC24 unless it proves to be a statistical outlier. (The January Septem- 
ber 2012 average of SSN is 59.5. The actual yearly value for SSN in 2012 turned out to be 57.5; the 
values in 2013 for January-May average about 62.) 

The inferred fit for ASAT (Rmax year) versus ASAT (Rmin year), having r- -0.295, is not 
statistically important ( cl < 90%); however, the inferred fit for GLOTI (Rmax year) versus GLOTI 
(Rmin year) is statistically important, having r- 0.798, se = (). 1 5 °C and cl> 99.5%. Applying the 
GLOTI (Rmin year) = 0.44 °C, one estimates that GLOTI (Rmax year) = 0.42 ±0.1 5 °C, or GLOTI 
(Rmax year) <0.57 °C, for SC24. For comparison, SC23 had GLOTI (Rmax year) = 0.36 °C. Inter- 
estingly, both ASAT and GLOTI values for the Rmin year are the warmest values on record for 
an Rmin year, being 9.75 °C and 0.44 °C, respectively. Perhaps this is a strong indication that both 
ASAT and GLOTI in the Rmax year for SC24 will exceed respective record values, being 10.07 °C in 
SC22 for ASAT and 0.36 °C in SC23 for GLOTI. 

Figure 13 shows the cyclic secular trend over SC9-SC23 for (a) ASAT (Rmin year), (b) GLOTI 
(Rmin year), (c) ASAT (Rmax year), and (d) GLOTI (Rmax year). While the inferred regressions for 
ASAT are not statistically important, those based on Fisher’s exact test do suggest a real secular rise 
in ASAT (Rmin year) and ASAT (Rmax year) over time. Hence, values of ASAT (Rmin year) and 
ASAT (Rmax year) appear to be slowly rising with the passage of time such that SC24 and later cycles 
more often than not should have ASAT (Rmin year) and ASAT (Rmax year) values warmer than 
their median values. In fact, since SC 17, six of eight SCs have had ASAT (Rmin year) values warmer 
than 9.23 °C, and since SC 16, seven of eight cycles have had ASAT (Rmax year) values warmer than 
9.2 °C. Based on the inferred regressions, one expects ASAT (Rmax year) = 9.67 ± 0.59 °C (or ASAT 
(Rmax year) <10.26 °C) for SC24 and ASAT (Rmin year) = 9.37 ±0.45 °C (or ASAT (Rmin year) 
<9.82 °C) for SC25. 
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Figure 13. Cyclic secular trend over SC9-SC23 of (a) ASAT (Rmin year), (b) GLOTI 
(Rmin year), (c) ASAT (Rmax year), and (d) GLOTI (Rmax year). 


For GLOTI, both the inferred regressions and Fisher’s exact tests reveal GLOTI (Rmin year) 
and GLOTI (Rmax year) to be statistically important, suggesting warmer values with the passage 
of time. For SC24, one estimates GLOTI (Rmax year) = 0.34 ±0.1 °C, (or GLOTI (Rmax year) 
<0.44 °C) and for SC25, one estimates GLOTI (Rmin year) = 0.32 ± 0. 16 °C (or GLOTI (Rmin year) 
<0.48 °C). 


24 


Figure 14 shows the scatter plots of (a) AS AT (Rmin year) and (b) GLOTI (Rmin year) 
versus Rmin. Neither parameter is found to display a statistically meaningful association. So, while 
the cyclic secular trends (see fig. 13) provide some indication as to the values of AS AT (Rmin year) 
and GLOTI (Rmin year), the actual value of the Rmin itself provides no clear indication as to the 
expected values of ASAT or GLOTI for the Rmin year. 



Figure 14. Scatter plots of (a) ASAT (Rmin year) and (b) GLOTI (Rmin year) versus Rmin. 


Figure 15 depicts the scatter plots of (a) ASAT (Rmax year) and (b) GLOTI (Rmax year) 
versus Rmax. While the inferred linear correlation between ASAT (Rmax year) and Rmax is only 
of marginal statistical importance (r = 0.36, cl > 90%), the inferred linear regression between GLOTI 
(Rmax year) and Rmax is highly statistically important (r = 0.67, cl > 98%). Presuming Rmax = 70 for 
SC24, one infers ASAT (Rmax year) = 8.93 ± 0.46 °C, or ASAT (Rmax year) <9.39 °C, and GLOTI 
(Rmax year) = -0.23 ±0.19 °C, or GLOTI (Rmax) <-0.04 °C. It is noteworthy that all previous 
cycles having Rmax <1 10.2 (the median Rmax based on SC12-SC23) have had GLOTI (Rmax year) 
<-0.03 °C. Hence, one strongly expects GLOTI (Rmax year) to be <-0.04 °C based on the expected 
Rmax = 70 (or smaller) for SC24, unless, of course, SC24 proves to be a statistical outlier, which now 
seems highly likely (based on GLOTI values for 2012 and expected values for 2013). 
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Figure 15. Scatter plots of (a) ASAT (Rmax year) and (b) GLOTI (Rmax year) versus Rmax. 


Regarding the inferred correlation between GLOTI (Rmax year) and Rmax, even though the 
fit is considered to be highly statistically important, 110 it is noteworthy that, relative to its regres- 
sion value (-0.04 °C based on Rmax= 119.6), SC23’s GLOTI (Rmax year) is about 2.1 se warmer 
(=0.36 °C, the warmest observed value in the modern era); therefore, if SC24’s GLOTI (Rmax year) 
behavior proves similar to that of SC23’s, then one cannot discount the possibility that GLOTI 
(Rmax year) for SC24 will be a positive value, perhaps considerably warmer than its inferred value 
of GLOTI (Rmax year) =-0.23 ±0.19 °C, using Rmax = 70. 


2.4 Behavioral Aspects of ASAT and GLOTI Relative to the AMO Index and the Trend 
in the Atmospheric Concentration of Carbon Dioxide at Mauna Loa, Hawaii 

As noted in section 2.3, factors other than solar forcing now appear to be the stronger fac- 
tors influencing the Earth’s climate system, at least over the past few decades. Examined here are the 
relations between the ASAT and GLOTI values against the AMO index values and the trend in the 
atmospheric concentration of carbon dioxide (C0 2 ) values as measured at Mauna Loa, Hawaii, for 
the interval 1950-2011, an interval that includes the current period when solar forcing of climate 
appears to have weakened. 


26 


The AMO is a temperature oscillation in the detrended sea surface temperature of the North 
Atlantic Ocean (0-70° N.) having a period of about 65-70 years that fluctuates between warm (posi- 
tive) and cool (negative) phases. 111 ’ 112 It is presently in its warm phase, having begun about 1995 
based on 10-yma values, with the minimum 10-yma value of the AMO index prior to the present 
warming phase occurring about 1975 (i.e., although the warm phase began about 1995, the trend 
in the warming began about 1975). The warmest 10-yma value of the AMO index for the current 
warm phase presently measures about 0.21 °C and is about 0.06 °C warmer than the preceding peak 
10-yma value that was seen in the last warming period (about 1956). It is warmer by about 0. 1 1 °C 
than the peak 10-yma value that was seen about 1876. The current warming phase is expected to 
last at least another one to two decades. Values of the monthly AMO index are available online at 
<http://www.esrl.noaa.gov/psd/data/correlation/amon.us.long.data>. 

Carbon dioxide is the major constituent of the greenhouse gases in the Earth’s atmo- 
sphere, 113-1 16 although other anthropogenic gases contribute as well, such as methane and nitrous 
oxide, but to a lesser extent as compared to C0 2 . 117 Continuous measurement of the atmospheric 
concentration of C0 2 has been ongoing at Mauna Loa since 1958. 118-120 While the atmospheric 
concentration of C0 2 is found to vary seasonally over the year, expressed as an annual average, it 
has been increasing unabatedly year-to-year from a mean of 315.97 ppm in 1959 to 391.57 ppm in 
2011, yielding an average rate of increase of about 1.45 ppm per year or about 14.54 ppm per decade. 
Annual averages of the Mauna Loa C0 2 (MLC02) measurements, expressed in ppm, are available 
online for the interval 1959-201 1 at <ftp://ftp.cmdl.noaa.gov/ccg/co2/trends/co2_annmean_mlo.txt> 
and annual values of the global radiative forcing due to greenhouse gases are available online for 
the interval 1979-2011 at <http://www.esrl.noaa.gov/gmd/aggi/>. (MLC02 measured 393.82 ppm 
in 2012.) 

Figure 16 shows the annual (thin jagged line) and 10-yma (thick smoothed line) of (a) the 
AMO index and (b) the MLC02 values for the interval 1950-2011. Regarding the variation of 
figure 16(a), it was reflective of the cool phase (negative 10-yma values) during the interval 1965- 
1994 and the warm phase (positive 10-yma values) both prior to 1965 and after 1994. As noted pre- 
viously, the present warming, as deduced using the 10-yma values, measured at least 0.21 °C in 2002 
and 2006, up from the cooling, which measured -0.25 °C in 1974 and 1975. The previous warming 
peaked at 0. 15 °C in 1956. (The current warming is expected to continue for at least another 20 years, 
possibly exceeding 0.21 °C.) 
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Figure 16. Annual (thin jagged line) and 10-yma (thick smoothed line) of (a) the 
AMO index and (b) the MLC02 values for the interval 1950-2011. 


Regarding figure 16(b), visually, the annual and 10-yma values are virtually identical. The 
10-yma value of MLC02 measures about 320 ppm in 1964 and about 382 ppm in 2006, an increase 
of about 19% in 42 years. Hence, based on the trend, atmospheric concentration of C0 2 at Mauna 
Loa is growing at an average rate of about 1.48 ppm per year, or 14.8 ppm per decade, with no evi- 
dence suggesting a leveling off. In fact, placement of a straight-edge along the curve shows that the 
MLC02 values actually are increasing faster than at a linear rate, which should be of great concern, 
especially because of the predicted consequences directly attributable to increased anthropogenic 
gases in the Earth’s atmosphere (and the associated warming of the atmosphere and oceans). 121-128 

Figure 17 displays the scatter plots of (a) AS AT and (b) GLOTI against the AMO index for 
the interval 1950-2006 based on 10-yma values. Two fits are shown in each subpanel with selected 
years identified (1950, 1975, 1976, and 2006). The lower fit is based on the interval 1950-1975, where 
the year 1975 represents the coolest extreme in the AMO index, and the upper fit is based on the 
interval 1976-2006, which is the current interval of continuous warming from its coolest AMO index 
value in 1975. Concerning figure 17(a), both fits are found to be statistically important, with the 
lower fit having r = 0.85 and cl> 99.9% and the upper fit having r = 0.97 and cl» 99.9%. Concerning 
figure 17(b), only the upper fit is found to be statistically important, having r = 0.99 and c/» 99.9%. 
Apparently, the ASAT and GLOTI values are now more strongly coupled to the AMO than they 
were prior to 1976. (The 10-yma values for the year 2007 for ASAT and GLOTI are 10.01 °C and 
0.59 °C, respectively. The 10-yma value of AMO for 2007 is 0.206 °C.) 
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Figure 17. Scatter plots of (a) AS AT and (b) GLOTI against the AMO index for the interval 
1950-2006, based on 10-yma values. 


Figure 18 shows the scatter plots of (a) AS AT and (b) GLOTI versus the MLC02 values for 
the interval 1964-2006 based on 10-yma values. Concerning figure 18(a), using all available values 
(denoted >’ 1964 - 2006 )’ the fit, while appearing to be statistically important, is found not to be statisti- 
cally important (r = 0.91, cl< 90%). However, when the comparison is based on the shortened inter- 
val 1982-2003 (denoted T 1982 - 2003 )’ the association between ASAT and MLC02 is indeed found to 
be statistically important (r = 0.99, cl >99. 9%). Concerning figure 18(b), it is inferred to be highly 
statistically important (r = 0.99, c/» 99.9%). (The 10-yma of MLC02 for 2012 is 383.67 ppm.) 



Figure 18. Scatter plots of (a) ASAT and (b) GLOTI versus the MLC02 values for the interval 
1964-2006, based on 10-yma values. 
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3. DISCUSSION AND CONCLUSIONS 


This TP has shown (see fig. 1) that the annual temperature as measured at the Armagh Obser- 
vatory and the temperature anomaly as described by the GLOTI have both revealed a warming 
planet Earth over the past century or more, with temperatures increasing at the average rate of 
0.14 °C per decade based on the ASAT and 0.09 °C per decade based on the GLOTI. 81 ’ 82 ’ 129 ’ 130 
Although the two parameters both show general long term secular increases over time, their individ- 
ual behaviors sometimes have differed. For example, the rise in temperature between about 1980 and 
2007 at Armagh appears to have been much steeper than was seen for the GLOTI. Likewise, while 
temperature was trending downwards between about 1945 and 1980 at Armagh, it is perceived to have 
been relatively flat as measured by the GLOTI during the same interval. Additionally, because the 
ASAT record begins in 1844 and the GLOTI record does not begin until 1880, the ASAT values have 
an advantage over the GLOTI values for determining temperatures prior to 1880. During this earlier 
epoch, the temperature at Armagh actually is found to have been warmer than its long term mean + 1 
sd of 9.25 ±0.55 °C and clearly was trending downwards from about 1845 to 1880. While true, the 
yearly values of the ASAT now being experienced are among the warmest ever recorded at Armagh, 
with the trend line now being warmer than was seen previously in 1849 by more than 0.5 °C, surpass- 
ing the 10-yma values of 1849 about the year 1990. 

Also shown in this investigation is that, while the behavior of SSN, ASAT, and GLOTI are 
strikingly similar, especially prior to SC21, the trends in temperature and SSN now no longer appear 
to be closely coupled. Hence, one surmises that other factors must now be dominating the Earth’s cli- 
mate system, at least over the past three SCs (i.e. , since SC21). 131 ’ 132 The trend in the GLOTI values 
remains upwards in 2006 (the last available 10-yma value), while the trend in the SSN values quite 
noticeably is downwards. The trend in the ASAT values appears to have peaked locally in 2002/2003, 
cooling slightly by 0.13 °C between 2003 and 2006. While true, it is anticipated that the trend in the 
ASAT values soon will flatten and then again continue to rise thereafter, especially given that the 
AMO remains in the warm phase and C0 2 continues to build up in the Earth’s atmosphere. 

Clearly, the ASAT and GLOTI values track each other very well, with 70% of the yearly val- 
ues obeying the paradigm of warmer (cooler) temperatures at Armagh being indicative of a warmer 
(cooler) temperature anomaly as measured by the GLOTI. However, a strange behavior (see fig. 4) 
has been found involving the 10-yma values of the two parameters. Namely, two separate regres- 
sions are inferred to link the GLOTI and ASAT values. The first regression spans 1885-1945, hav- 
ing r = 0.85, r 2 = 0.72, and se = 0.06 °C. The second regression spans the current interval 1982-2006, 
having r = 0.96, r 2 = 0.93, and se- 0.04 °C. The interval between 1945 and 1982 appears to be one of 
transition with 10-yma values of ASAT displaying cooling and 10-yma values of GLOTI displaying 
a near flat appearance. The two regressions are offset by about 0.3 °C, with the current values being 
warmer. Hence, while a 10-yma value of ASAT equal to 9.5 °C as measured in 1942 yields a 10-yma 
GLOTI of about -0.05 °C, such a value now yields a 10-yma value of GLOTI equal to about 0.27 °C, 
a difference of 0.32 °C. It is unclear as to what may have caused the relationship between the 10-yma 
values of GLOTI and ASAT to change and be displaced upwards by about 0.32 °C. 
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This TP also has shown that the cyclic variation of the <ASAT> and <GLOTI> has been 
decidedly upwards towards warmer temperatures since SC 12, with no evidence of flattening or cool- 
ing except for a brief interval between SC18 and SC21 (see fig. 7). While cyclic averages of <SSN> 
clearly are now tracking downwards (since SC21), cyclic averages of <ASAT> and <GLOTI> remain 
tracking upwards towards warmer temperatures. Hence, solar forcing of Earth’s climate system, 
which should result in cooling if it remained the main forcing factor in climatic change, now seems to 
have been overwhelmed by other factors, factors that have resulted in a continued warming as shown 
by the ASAT and GLOTI values, even as solar activity has declined. Clearly, SC23 was warmer than 
SC22, which was warmer than SC21. The inference then is that the temperature at Armagh possibly 
could be warmer in SC24 as compared to what it was in SC23 and that the <GLOTI> values possibly 
could be higher as well. 

Figure 19 replots figures 7(e) and (f) as linear fits of (a) <ASAT> and (b) <GLOTI> ver- 
sus SC for SC12-SC23. Both fits are inferred to be highly statistically important. Using the pre- 
ferred fits (y) based on SC12-SC23, one estimates SC24 to likely have <ASAT> = 9.82 ± 0.24 °C and 
<GLOTI> = 0.37 ±0.09 °C (the ±1 se prediction intervals), inferring only about a 16% chance that 
<ASAT> and <GLOTI> for SC24 will be >10.06 °C and >0.46 °C, respectively. Hence, one easily 
could conclude that it appears unlikely that SC24 will have <ASAT> and <GLOTI> values warmer 
than was seen in SC23 (i.e. , 10.08 °C and 0.49 °C, respectively), as suggested by Solheim et al. 61 
(In figure 19, the error bars represent the mean ± 1 sd interval for each cycle.) 



Figure 19. Finear fits of (a) <ASAT> and (b) <GFOTI> versus SC for SC12-SC23. 
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However, alternate regressions (y') plotted in figure 19 based only on the most recent SCs 
(SC20-SC23) counter this conclusion, where these more recent cycles suggest a faster warming from 
one cycle to the next than was seen in the earlier cycles. While the inferred regression for <ASAT> 
does not appear to be statistically important, the one for <GLOTI>, in contrast, is found to be highly 
statistically important (r = 0.99, cl >99%). Hence, presuming the validity of the alternate regression 
fit, <GLOTI> is inferred to be about 0.64 ±0.03 °C for SC24, thus suggesting that SC24 undoubt- 
edly will be warmer than SC23, with only about a 16% chance of being cooler than 0.61 °C, a value 
at least 0. 12 °C warmer than was seen in SC23 (0.49 °C). 

It should be noted that in this analysis, SC23 (1996-2007) stands out as having <ASAT> and 
<GLOTI> values considerably warmer (the warmest of any modern era SC on record) than would 
have been expected based on its observed Rmax or <SSN> (see fig. 8). Likewise, most unusual are 
the values of the AS AT and GLOTI for the year 2010, a year when the AS AT had the lowest tem- 
perature (8.75 °C) since 1986, yet the GLOTI had the highest temperature anomaly (0.63 °C). 

For SC24’s Rmax year (believed to be either the year 2012 or 2013), the inferred regression of 
GLOTI (Rmax year) versus SC suggests a value of about 0.34 ±0.1 °C (the ±1 se prediction inter- 
val) for GLOTI (see fig. 13). Based on its expected Rmax value of about 70, GLOTI (Rmax year) 
is expected to be only about -0.23 ±0.19 °C (see fig. 15). Presently, the GLOTI values for 2012 are 
averaging about 0.5 °C (based on the January-September average). Hence, if the year 2012 proves 
to be the Rmax year, it will be the warmest such year on record, exceeding that measured for SC23 
(0.46 °C in the year 2000). (For GLOTI, the year 2012 measured 0.56 °C.) 

Comparisons of 10-yma values of the AS AT and GLOTI against the AMO index show that 
the current interval 1976-2006 (SC2 1-present) has inferred regressions that are highly statistically 
important (r = 0.97 and cl» 99.9% for ASAT and r = 0.99 and cl»99.9% for GLOTI), suggesting that 
the current interval is one where the ASAT and GLOTI values are highly coupled with the AMO. 
Ten-yma values of the AMO index values exceeding 0.2 °C are suggestive that 10-yma values of 
ASAT >10 °C and GLOTI values >0.5 °C. (The AMO index currently reflects a continuing warming 
phase, one that is expected to persist for at least 10-20 years or more. Furthermore, the 10-yma value 
of AMO exceeded 0.2 °C in 2006 and 2007.) 

Comparisons of 10-yma values of the ASAT and GLOTI against the 10-yma values in 
MLC02 for the interval 1964-2006 show that, while the inferred regression of ASAT versus MLC02 
does not appear to be statistically important, the one for the GLOTI versus MLC02 is inferred to 
be highly statistically important (r = 0.99 and cl» 99.9%). Because the trend in MLC02 values clearly 
is upwards, increasing at the average rate of about 1.48 ppm per year based on the 10-yma values, 
atmospheric concentration of C0 2 at Mauna Loa probably will exceed 400 ppm by the year 2015. 
A 10-yma value of MLC02 equal to 400 ppm suggests a 10-yma value of GLOTI equal to about 
0.72 ±0.06 °C (see fig. 18). 

In conclusion, the trend in temperature at Armagh now appears to be strongly coupled with 
the trend in the AMO. Likewise, the GLOTI now appears to be strongly coupled to both the trends in 
the AMO and the level of atmospheric concentration of C0 2 as measured at Mauna Loa. Since the 
trend in AMO remains reflective of the warm phase, which is expected to persist for the next decade 
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or longer, the trend in yearly temperatures at Armagh is expected to continue to be warmer than the 
long term average (9.25 ±0.55 °C) and may even rise to new record highs (i.e., annually >10.6 °C, 
10-yma trend value >10.13 °C). Similarly, the trend in the GLOTI values is expected to continue to 
rise, especially because of the continuing rise in the trend of the anthropogenic gases in the Earth’s 
atmosphere. As to whether the < ASAT> and <GLOTI> values will be warmer in SC24 as compared 
to that seen in SC23 (i.e., 10.08 °C for <ASAT> and 0.49 °C for <GLOTI>), presently, one cannot 
be completely certain. Based on the first 4 years of SC24 (2008-201 1), the ASAT is averaging 9.65 
±0.65 °C in SC24, as compared to 9.96 ±0.5 °C in the first 4 years of SC23 (1996-1999), while the 
GLOTI is averaging 0.54 ±0.08 °C in the first 4 years of SC24, as compared to 0.42 ±0.1 3 °C in the 
first 4 years of SC23. 
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